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Abstract — Optimal  waveforms  for  minimum  mean 
square  error  range  profile  estimation  are  investigated. 
An  idealized  measurement  and  waveform  adaptation 
process  is  developed  that  yields  optimal  scene  and 
range  specific  waveforms.  This  process  is  idealized  in 
that  during  each  cycle  of  the  process,  a  large  number 
of  dwells  are  required.  As  part  of  our  method,  a 
modified  version  of  the  Adaptive  Pulse  Compression 
(APC)  estimation  method  is  used  to  estimate  the 
range  profile  after  each  dwell  cycle.  The  proposed 
method  is  analogous  to  the  APC  method  in  that  it 
yields  a  set  of  range  specific  optimal  waveforms,  while 
the  APC  method  yields  a  set  of  range  specific  optimal 
pulse  compression  filters.  In  certain  scenarios,  the 
measurement  and  waveform  adaptation  process  yields 
range  profile  estimates  that  are  significantly  better 
than  those  derived  by  the  APC  method  alone. 

1.  Introduction 

In  modern  radar  and  other  active  sensors,  a  signal  is 
transmitted,  directed  into  the  associated  medium  such 
as  free  space  in  radar,  and  reflected  back  to  a  receiver 
by  any  scatterers.  The  delays  associated  with  reflected 
signals  arriving  at  the  receiver  provide  a  spatial  profile 
associated  with  the  respective  ranges  to  the  scatterers. 

In  radar,  the  high-power  narrow-pulse  transmissions 
required  to  achieve  high  range  resolutions  are  often 
limited  by  practical  power  constraints.  To  overcome  such 
limits,  transmitted  radar  signals  often  take  the  form  of 
a  phase-  or  frequency-  modulated  rectangular  pulse  to 
obtain  high  spatial  resolution  inversely  proportional  to 
the  modulation  bandwidth  in  a  technique  known  as  pulse 
compression  [4].  The  received  signal  consists  of  delayed, 
attenuated  versions  of  the  transmitted  waveform.  A  Alter 
matched  to  the  transmitted  waveform  extracts  the  high- 
resolution  spatial  profile  from  the  received  signal. 

For  a  solitary  point  scatterer  in  the  presence  of  white 
Gaussian  noise,  the  matched  Alter  (i.e.,  matched  to  the 
transmitted  signal  and  system  noise)  maximizes  output 
signal-to-noise  ratio  (SNR)  and  the  detectability  of  scat¬ 
terers  in  noise.  The  matched  Alter  correlates  the  delayed, 
attenuated  versions  of  the  transmitted  waveform  as  they 
arrive  at  the  receiver  and  provides  the  range  and  relative 
complex  amplitudes  of  scatterers.  The  correlation  of  the 
matched  Alter  with  the  waveform  is  the  autocorrelation 
of  the  waveform.  However,  an  inherent  problem  is  the 
masking  of  small  targets  by  large  nearby  targets  due  to 


range  sidelobes  that  are  present  when  using  the  standard 
matched  Altering  method. 

Blunt  and  Gerlach  [1]  developed  a  method  for  range 
proflle  estimation.  Adaptive  Pulse  Gompression  (APG), 
in  which  a  Alter  adapted  to  the  scene  is  derived  for 
each  range  bin.  Each  Alter  is  optimally  adapted  to  its 
respective  range  cell  based  on  the  minimum  mean-square 
error  (MMSE)  criterion  [2].  The  adaptation  process  for 
deriving  optimal  Alters  is  iterative  in  the  sense  that  a 
numerical  solution  method  is  iterative.  However,  it  only 
requires  the  result  of  a  single  radar  dwell  as  input. 

The  adapting  of  the  signal  processing  with  respect 
to  the  range  bin  allows  the  sidelobes  of  targets  to  be 
greatly  reduced,  thereby  mitigating  to  a  large  extent 
the  aforementioned  problem  of  the  masking  of  small 
targets  by  nearby  larger  targets.  Speciflcally,  APG  has 
been  demonstrated  to  mitigate  autocorrelation  sidelobes 
to  the  level  of  the  noise  and  accurately  estimate  small 
targets  in  the  presence  of  several  larger  targets.  Gompar- 
isons  of  this  Adaptive  Pulse  Gompression  (APG)  method 
and  the  standard  matched  Alter  estimation  showed  it  to 
be  superior  over  a  variety  of  stressing  scenarios. 

Despite  the  fact  that  APG  derives  optimal  Alters  for 
each  range  cell,  for  radar  scenes  in  which  there  are  groups 
of  closely  spaced  targets,  even  after  APG  processing  there 
can  be  signiflcant  sidelobes  near  the  group  of  targets 
that  could  potentially  mask  targets  close  to  the  cluster  of 
targets.  Such  a  radar  scene  could  for  example  represent 
a  littoral  environment.  In  this  paper,  the  APG  method 
is  extended  by  deriving  optimal  scene  and  range-speciflc 
waveforms.  These  optimal  waveforms,  in  conjunction 
with  optimal  Alter  weights,  have  been  found  in  simula¬ 
tion  to  further  enhance  the  estimation  performance  for 
certain  scenarios  with  clusters  of  targets. 

The  framework  of  APG  is  utilized  here.  The  signal 
model,  measurement  model,  and  performance  criterion 
are  the  same  as  those  utilized  in  APG.  Also  as  in  APG, 
the  waveforms  are  restricted  to  being  coded  waveforms 
with  a  Axed  number  of  chips.  Finally,  the  derivation  for 
optimal  waveforms  is  analogous  to  the  APG  derivation 
for  optimal  weights.  Accordingly,  the  space  of  admissible 
waveforms  over  which  optimization  is  done  is  such  that 
for  any  particular  waveform,  the  chips  of  the  waveform 
do  not  all  necessarily  have  the  same  amplitude.  While 
the  derivations  are  mathematically  parallel  to  each  other. 
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the  measurement  and  adaptation  process  required  differs 
greatly. 

As  mentioned  previously,  in  APC,  the  adaptation  pro¬ 
cess,  while  iterative,  only  requires  as  input  the  result  of  a 
single  dwell.  This  is  because  of  the  effect  of  changing  a  fil¬ 
ter  can  be  observed  by  simply  processing  the  radar  return 
again  using  the  new  filter.  In  contrast,  the  adaptation 
process  for  waveforms  requires  multiple  dwells  since  the 
effect  of  the  change  in  a  waveform  can  only  be  observed 
by  transmitting  another  dwell. 

We  have  posited  a  joint  measurement  and  adaptation 
process  that  assumes  range  cell  specific  waveforms  are 
repeatedly  tested  against  their  respective  range  cells  by 
transmitting  them  during  the  dwell  cycle  of  the  pro¬ 
cess.  That  is,  if  there  are  L  range  cells,  then  multiple 
sets  of  L  waveforms  are  transmitted  during  each  cy¬ 
cle  of  the  process.  This  process  is  probably  unrealistic 
except  for  certain  specific  situations,  and  as  such  is 
idealized.  However,  the  waveforms  derived  using  it  are 
best-case  waveforms,  and  thus  given  an  indication  of 
the  upper  bound  on  the  estimation  performance  that 
can  be  achieved  using  joint  waveform/filter  adaptation. 
Additionally,  this  process  could  potentially  be  used  as 
the  basis  of  measurement  processes  for  specific  scenarios. 
For  example,  a  modified  process  could  be  useful  for  a 
target  tracking  scenario  where  only  a  few  range  bins  are 
of  interest  as  in  [3] .  Another  example  could  be  a  scenario 
where  there  is  existing  knowledge  of  the  environment. 

We  point  out  that  a  modified  version  of  APC  is  used 
to  estimate  the  range  profile  as  part  of  the  measurement 
and  waveform  adaptation  process.  Specifically,  the  APC 
algorithm  is  used  to  estimate  the  range  profile  using  the 
return  from  a  different  waveform  for  each  range  cell. 
The  measurement  and  adaptation  process  was  simulated 
for  a  few  scenarios,  and  was  found  to  yield  estimation 
performance  superior  to  that  resulting  from  APC  alone. 

An  agile  waveform  method  has  been  proposed  that 
derives  optimal  waveforms  in  heavy  sea  clutter  in  [3]. 
In  [3],  sidelobes  of  nearby  scatterers  are  explicitly  mini¬ 
mized  at  a  particular  location  of  interest,  in  the  context 
of  target  tracking.  Here,  the  sidelobe  minimization  is 
a  by-product  of  the  minimization  of  the  mean-square 
error.  Furthermore,  in  [3],  the  waveforms  are  restricted  to 
constant  amplitude  waveforms,  while  here,  less  practical 
amplitude  modulated  waveforms  are  allowed. 

II.  Waveform  Adaptation 

Our  notation  is  the  same  as  that  of  [I],  with  minor 
exceptions  as  noted.  For  a  more  detailed  exposition, 
please  refer  to  [1].  The  signal  model  assumed  here  is 
the  same  one  that  was  defined  in  [1],  except  that  we 
assume  a  range-specific  waveform  s{l)  is  used  for  each 
range  1  <  I  <  L.  The  range  profile  to  be  estimated  is 
X  =  [  a;(l)  x{2)  . . .  x{L)  ]  where  E[x{l)x{j)]  =  0 

ior  I  ^  j ,  1  <  I  <  L,  1  <  j  <  L.  The  observation  resulting 
from  using  waveform  s(l)  is  assumed  to  be  given  by 

y(0  =  A^(0s(/)+v(0  (1) 


where  v(^)  is  the  measurement  noise  and  A{1)  is  an 
N  X  N  matrix  where  the  jth  row,  1  <  j  <  fV,  is  given  by 
[x{l-j  +  l)  x{l-j  +  2)  ...  x{l-j  +  N)]. 

Denote  the  range  profile  estimate  by  x  = 

[  a:(l)  x(2)  . . .  x{L)  ] .  The  estimate  x{l)  of 
x{l)  is  of  the  form 

x{l)  =  w^{l)y{l),  (2) 

and  consequently  the  estimation  objective  function  to  be 
minimized  via  choice  of  waveforms  and  adaptive  weights 
is,  as  defined  in  [1], 

J(0  =  A[|x(0-w«(0y(0P] 

=  E[{x{l)-w{l)^y){x{ir-y^wm 
=  E[{x{l)x{l)*  -  w{l)^y{l)x*{l)  -  x(0y^(/)w(/) 

+w^{l)y{l)y^{l)w{l)] 

=  E[x{l)x*{l) -w^{l)A^{l)s{l)x*{l) 

— w^(/)v(?)x*(/)  —  x{l)s^  {l)A*{l)w{l) 
—x{l)v{l)^'w{l)  +  {l)A^s{l)s^  {l)A*'w{l) 

{l)A^  {l)s{l)v^  {l)w{l) 

{l)v{l)s^  {l)A*  {l)w{l) 

+'w^  {l)w{l)] 

Setting  d.J{l)/ds^{l)  =  0  yields  that  for  given  w(l),  the 
waveform  s(/)  that  minimizes  J{1)  satisfies 

A[A*(Z)w(Z)w(Z)^A^(0]s(0  =  A[a;(0A*(0]w(0 

Using  the  fact  that  the  range  cells  are  uncorrelated  with 
each  other  yields 

E[A*{l)w{l)'w^  {l)A'^{l)]s{l)  =  p{l)w{l) 

where  p{l)  =  E[x*  {l)x{l)].  By  using  this  fact  again,  this 
equation  becomes 

N-l 

(  H  P{l-n)wn{l)w^{l))s{l)  =  p{l)w{l)  (3) 

n=—N-\-l 

where  w„(?)  contains  the  elements  of  the  vector  of 
weights  w(^)  =  [  wo{l)  wi{l)  ...  wn-i{1)  ] 

shifted  by  n  samples  and  zero-fillech  e.g. 
W2(0  =  [00  Wo{l)  ...  WN-sil)]  and 

w_2(0  =  [  1x2(1)  ...  wn-i(1)  00].  When 

(I2n=-N+iP(^  ~  exists,  the  optimal 

waveform  for  range  cell  I  can  be  expressed  as 

N-l 

®(o  =  (  p(^-^)w«(0wf(0)”V(0w(0  (4) 

As  was  shown  in  [1],  the  filter  w(^)  that  minimizes  J(l) 
for  a  given  waveform  s(l)  is  given  by 

N-l 

w(0  =  p(/)(R+  ^  p{l  +  n)s^{l)s^{l))-h{l)  (5) 

where  R  =  E[v{1)n^ (1)].  We  note  that  p{l)  is  estimated 
simply  by  p(l)  =  x*(l)x(l). 
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It  is  desired  to  find  the  waveform  and  filter/ weights 
that  minimize  J{1)  for  each  range  cell.  However,  as  men¬ 
tioned  previously,  the  effect  of  updating  a  waveform  can 
only  be  observed  by  actually  transmitting  the  waveform. 
The  measurement  and  adaptation  process  for  deriving 
optimal  waveforms  is  thus  measurement  intensive  in 
contrast  with  APC,  which,  as  it  uses  a  common  waveform 
for  all  range  cells,  only  requires  the  result  of  a  single 
measurement/radar  dwell. 

In  our  measurement  and  adaptation  process,  two  cy¬ 
cles,  a  waveform  update  cycle  (WC)  and  a  filter  update 
cycle  (FC),  are  repeatedly  performed  one  after  the  other. 
The  filter  remains  fixed  during  the  waveform  update 
cycle  (WC),  and  the  waveform  remain  fixed  during  the 
filter  update  cycle  (FC). 

The  waveform  update  cycle  (WC)  consists  of  several 
sets  of  dwells.  During  each  set  of  dwells,  each  range 
cell  is  interrogated  by  a  waveform  specific  to  that  range 
cell,  resulting  in  the  observation  (1).  After  all  range- 
specific  waveforms  in  a  set  have  been  transmitted,  the 
range  profile  estimate  and  waveforms  are  updated,  and 
then  this  process  repeats,  i.e.,  another  set  of  dwells  is 
transmitted.  The  range  profile  estimate  is  updated  by 
applying  the  current  filter  w(/)  to  the  observation  y(0 
according  to  equation  (2).  The  waveforms  are  updated 
by  solving  (3)  (via  matrix  inversion)  using  the  updated 
range  profile  estimate  and  the  current  range  specific 
weights. 

During  the  filter  update  cycle  (FC),  the  weights  are 
updated  using  APC,  but  with  a  waveform  specific  to  each 
range  cell,  instead  of  a  common  waveform  for  all  range 
cells.  Specifically,  in  each  cycle  of  the  modified  APC 
algorithm,  the  range  profile  estimate  is  first  updated 
according  to  (2)  for  each  1  <  I  <  L,  and  then  the  weights 
•w{l)  are  updated  according  to  (5)  for  all  1  <  ?  <  L. 

Cycles  of  waveform  updates  (WC)  followed  by  filter 
updates  (FC)  are  repeated  until  the  range  profile  esti¬ 
mate  converges,  which  in  practice  is  after  only  a  few  fixed 
number  of  cycles.  Of  course,  each  waveform  update  cycle 
consists  of  a  large  number  of  measurements/radar  dwells. 
As  in  the  APC  method,  the  estimates  x{l)  calculated  in 
either  a  waveform  update  cycle  or  a  filter  update  cycle  are 
bounded  below  in  magnitude  to  prevent  ill-conditioning 
in  the  matrix  inversion. 

III.  Simulation  Results 

In  this  section,  the  adaptive  measurement  process 
described  in  the  previous  section  is  simulated.  Two  sce¬ 
narios  are  simulated.  In  both  scenarios,  the  true  radar 
scene  is  combined  with  clutter  as  well  as  measurement 
noise  for  each  dwell  measurement.  We  note  that  our 
simulation  of  clutter  is  such  that  the  clutter  changes 
from  measurement  (dwell)  to  measurement,  as  the  large 
number  of  measurements  would  make  it  unrealistic  to 
assume  fixed  clutter  over  an  entire  set  of  measurements. 
Specifically,  random  clutter  is  generated  anew  for  each 
measurement  of  the  waveform  cycle  described  above.  The 


clutter  and  measurement  noise  levels  are  both  set  65dB 
down  from  the  peak  of  a  unity  amplitude  target.  For 
all  scenarios,  the  measurement  and  adaptation  process 
initially  uses  a  32  bit  P4  code  for  the  waveform  and 
matched  filter.  In  all  scenarios,  the  phases  of  the  targets 
are  assumed  to  be  random  and  uniformly  distributed. 

In  one  scenario,  the  radar  scene  assumed  for  the  mea¬ 
surement  and  adaptation  process  consists  of  a  “pedestal” 
of  targets  of  the  same  magnitude,  with  random  phases. 
For  this  scenario,  the  robustness  of  the  optimal  wave¬ 
forms  is  investigated  by  applying  the  waveforms  to  a 
different  scene  in  which  the  pedestal  alternates  between 
target  and  no  target  from  range  cell  to  range  cell. 

Fig.  1  shows  the  range  profile  estimate  that  results 
from  the  APC  method  alone.  Fig.  2  shows  the  range 
profile  estimate  resulting  from  the  measurement  and 
waveform  update  process.  From  these  figures,  it  is  appar¬ 
ent  that  “sidelobes”  from  the  cluster  of  targets  is  about 
lOdB  higher  for  APC  alone  than  for  the  joint  waveform 
and  weight  adaptation  process. 

Fig.  3  shows  the  optimal  waveforms  for  all  range  cells. 
Near  the  cluster  of  targets,  the  waveforms  start  to  ap¬ 
proach  an  impulsive  waveform,  but  still  have  significant 
though  reduced  amplitudes  away  from  their  main  spikes. 
Fig.  4  shows  the  optimal  filters  for  all  range  cells.  From 
this  figure  it  can  be  seen  that  the  filter  weights  are 
high  near  the  edges  of  the  target  cluster,  where  the 
waveforms  are  quite  “spiky”  (see  Fig.  3).  Fig.  5  shows  the 
filter  responses  for  each  range-specific  optimal  waveform- 
filter  pair.  The  responses  are  very  much  like  the  desired 
“thumbtack”  response,  except  near  the  edges  of  the  tar¬ 
get  cluster.  In  these  regions,  a  filter  is  orthogonal,  via  the 
adaptation,  to  a  composite  of  the  shifts  of  the  waveform 
corresponding  to  the  surrounding  clutter  environment.  It 
is  evident  in  Fig.  5  that  the  adaptive  filters  have  shifted 
the  range  sidelobes  of  the  waveform/filter  pairs  according 
to  the  surrounding  environment.  The  optimization  of 
waveforms  eases  the  requirement  of  the  adaptive  filter 
weights  for  a  particular  range  cell  to  mitigate  interference 
from  the  surrounding  environment. 

Fig.  6  shows  the  range  profile  estimate  resulting  from 
using  the  waveforms  and  filters  derived  for  the  “solid” 
cluster  of  targets  for  a  cluster  of  targets  with  spaces  be¬ 
tween  the  targets.  In  this  case,  the  waveform-filter  pairs 
apparently  perform  well  away  from  the  target  cluster, 
where  the  radar  scene  is  largely  unchanged  from  the 
scene  used  for  optimization.  In  particular,  a  small  target 
has  been  added  to  the  right  of  the  cluster  of  targets, 
and  is  visible  even  though  it  was  not  included  in  the 
original  radar  scene  used  for  optimization.  However,  the 
pairs  do  not  perform  as  well  within  the  spaces  in  the 
target  cluster,  since  in  the  optimization  scene  the  spaces 
were  targets. 

Fig.  7  shows  the  range  profile  estimate  resulting  from 
utilizing  for  all  range  cells  the  waveform-filter  pair  that 
is  optimal  for  a  location  to  the  left  of  the  cluster  where 
an  additional  small  target  has  been  added.  Apparently, 
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this  waveform-filter  pair  is  not  particularly  good  for  a 
number  of  range  cells,  due  to  the  range-specific  nature 
of  the  pair.  Fig.  8  shows  the  range  profile  estimate  using 
the  waveform  that  is  optimal  for  the  small  target,  but 
allowing  the  filters  to  be  optimized  over  all  range  bins 
using  APC.  In  this  case,  the  performance  is  significantly 
improved  over  Fig.  7.  Fig.  9  shows  the  magnitude  of  the 
waveform  used  in  Fig.  7  -  8.  As  can  be  seen,  this  waveform 
on  one  hand  is  quite  similar  to  an  impulsive  waveform, 
but  on  the  other  hand  it  still  has  significant  amplitude 
away  from  its  peak. 

Another  scenario  simulated  involves  a  cluster  of  targets 
of  various  amplitudes,  again  with  random  phases.  Fig. 
10  shows  the  range  profile  estimate  resulting  from  APC 
alone.  Fig.  12  shows  the  range  profile  estimate  resulting 
from  the  joint  measurement  and  waveform-filter  adapta¬ 
tion  process.  Again,  the  joint  process  results  in  a  range 
profile  estimate  that  is  better  than  that  of  the  APC 
method  alone.  Fig.  11  shows  the  range  profile  estimate 
resulting  from  a  modified  measurement  process  in  which 
only  the  waveforms  are  optimized.  The  weights  used  are 
the  original  matched  filter  weights,  except  that  they  are 
normalized  with  respect  to  each  range  specific  waveform. 
In  this  case,  the  resulting  range  profile  estimate  appears 
to  have  about  the  same  accuracy  as  the  APC  estimate. 
However,  the  peak  signal  levels  are  inaccurate  due  to 
the  fact  that  the  waveforms  are,  after  each  iteration, 
normalized  to  have  a  maximum  unity  amplitude. 

IV.  Conclusion 

An  measurement  and  adaptation  process  that  jointly 
updates  range  specific  waveforms  and  filters  to  optimize 
them  with  respect  to  the  minimum  mean  square  error 
criterion  was  developed.  The  process  was  simulated  for  a 
few  scenarios,  and  was  found  to  yield  significant  improve¬ 
ment  over  APC  alone.  The  robustness  of  the  resulting 
optimal  waveform-filter  pairs  was  examined  by  testing 
them  against  scenarios  derived  by  altering  the  scenario 
for  which  the  pairs  were  optimized.  Although  idealized, 
the  process  developed  could  potentially  be  adapted  for 
specific  scenarios  of  interest,  such  as  littoral  scenarios 
with  pre-existing  environmental  knowledge,  or  target 
tracking,  in  which  only  a  very  limited  number  of  range 
bins  are  of  interest.  Finding  optimal  waveforms  within 
the  class  of  phase  coded,  constant  amplitude  waveforms, 
using  the  APC  framework,  is  a  topic  for  future  research. 
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APC  Range  Profile  Estimate 


Optimal  Filters 


Range  Profile  Estimate  with  Waveform  Adaptation  after  Sets  of  Dwells  with  APC 


Optimal  Waveforms 


Fig.  3.  Optimal  waveforms  for  pedestal  target  scene 


Responses  for  Range-Specific  Optimal  Filter-Waveform  Pairs 


Range  Bin 


Fig.  5. 
pairs 


Time  Epoch 

Filter  responses  for  optimal  waveform-filter 


Range  Profile  Estimate  Using  Sub-Optimal  Waveform-Filter  Pairs  with  Small  Adjacent  Target 
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Range  Profile  Estimate  using  Single  Waveform-Filter  Pair  with  Adjacent  Small  Target 


adjacent  small  target 


APC  Range  Profile  Estimate 


Range  Profile  Estimate  Using  Sub-Optimal  Waveform  with  APC  Processing 


dwells  (without  APC)  for  multiple  target  scene 


Waveform  tor  Range  Bin  Near  the  Pedestal 


Range  Profile  Estimate  with  Waveform  Adaptation  after  Sets  of  Dwells  with  APC 


Fig.  12.  Range  profile  estimate  resulting  from  mea¬ 
surement  process  consisting  of  cycles  of  sets  of  dwells 
followed  by  APC  for  scene  with  multiple  targets 
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